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Abstract The continuum one-photon annihilation at 7(2S) in ete~ experiment is 
studied. Such contributions to the measured final state wx? and m+ at u(25) mass 
are estimated by phenomenological models. It is found that these contributions must 
be taken into account in the determination of branching ratios of (2S) — wr’ and 
W(2S8) — nr, as well as other electromagnetic decay modes. The study reaches the 
conclusion that in order for BES to obtain the correct branching ratios on these decay 
modes, at least 10pb~! of data below the w(2.9)) peak is needed. 
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1 Introduction 


The study and understanding of the decay dynamics of the charmonium states J/w, w(2S) 
and possibly ~(3770), is one of the most important topics in Quantum Chromodynamics 
(QCD). There have been lots of activities in experimental part of this study in recent 
years since BES has collected the world largest (2S') sample during 1993-1994 and 
1995 running years at the w(2S) peak. Among the studies, the reconfirmation of the 
vector-pseudoscalar (VP) decay puzzle (also called “pz puzzle”) between J/w and w(2S) 
decays first observed by MARK II f] and the first observation of the ~(25) decay 
suppressions in vector-tensor (VT) modes ff] compared to naive perturbative QCD pre- 
dictions are of great interest. Since then, many theoretical efforts have been put to figure 
out the possible reasons why J/w and y(25) decays to some final states, especially for 
VP and VT modes significantly deviate from “15 % rule”, while many other decay modes 
agree with the predictions [f]. 
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To extend the experimental study to other decay modes of wW(2S), to further explore 
the studied modes to higher sensitivity, and to match the huge sample BES has collected 
at the J/w peak during 1999-2001 running years [6], BES has collected 14 M 7~(2S) data 
during 2001-2002 running year f]. This world largest sample makes it possible to measure 
the w(2S) decays to electromagnetic decay level, so that one can reach the sensitivity of 
separating the QED and QCD amplitudes in ~(25S') decays. 


As for the study of (3770), as an attempt to solve the “pa puzzle” in w(2S) decays, it 
has been proposed that the suppression of the y(25) decays to some specific final states 
is due to the mixing between S- and D-wave charmonium states, or the ¥(2S) is very 
close to DD threshold so that virtual process of (29) — DD cancels lots of (2.5) decay 
amplitudes to those suppressed channels [7]. Furthermore, the study of the old data may 
indicate 7)(3770) has a large fraction of decays into light hadrons besides DD modes, in 
contradiction to the picture that it decays to DD dominantly [B]. To confirm this, the 
best way is to measure the DD and total 7(3770) cross sections with high precision, so 
that the difference between the two gives the contribution of the non-DD decay rate, 
while the other possible way is to measure the exclusive non-DD decay modes. 


Now we know that J/w and w(2S) decay into hadrons, besides those with charmonium 
in final particles, through two interactions: the one-photon electromagnetic interaction 
and three-gluon strong interaction. The amplitudes of them, in general, may have a 
relative phase. This is also true for w(3770) decaying into light hadrons. 


In ete colliding beam experiment, J/w and w(2S) are produced by e*e7~ annihilation, 
there is inevitable the process 
ete — * — hadrons (1) 
produced simultaneously, which is indistinguishable from the hadron events from J/q or 
w(2S) decays. So in ete~ experiment, any final state may come from three processes: 
the charmonium three-gluon decays, the charmonium one-photon decays, and the virtual 
photon process without going through resonance (continuum process). The study of the 
charmonium decay dynamics is to determine the decay amplitudes of the three-gluon 
and the one-photon processes. In the case of the w(2S), the one-photon decays has a 
comparable cross section with respected to the continuum process. Moreover, for some 
strongly suppessed modes like VP and VT, the three-gluon process may also have a 
comparable cross section. In this case, one has to consider three amplitudes and the 
relative phases between any of the two amplitudes in the analysis of the experimental 
data. Fortunately, there are some decay modes of the charmonium where the strong 
decays are forbidden, only decays through one-photon annihilation is allowed, like (2S) 
— wr, which violates the isospin, and ~)(29) — ntr, which violates the G-parity. In 
such situation, only two amplitudes and one relative phase are present, this substantially 
simplifies the study, and makes the determination of the one-photon decay amplitude 
possible, provided the amplitude of the continuum process being known. 


The electromagnetic processes, such as ete~ — wr? and ete” — ata, are similar 
to the process ete” — yt” in the way that there are two Feynman diagrams : the 


continuum one-photon diagram and the (2S) diagram. Taking n*m" final state as an 
example, two diagrams are shown in Fig. [I]. 
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Figure 1: Two Feynman diagrams of ete~ — ntr” final state at w(2S). 


The observed experimental cross sections of these exclusive channels consist of three 
parts: the w(2S) resonance, the continuum and their interference. However, unlike the 
utu channel which the continuum amplitude can be calculated by QED, such terms for 
wr and ttr are to be determined by experiment, i.e., we need to measure the form 
factors of wr? and mtr. 


It is important to notice that the measured cross sections at the narrow resonances are 
sensitive to the experimental conditions, but the three parts in the cross section depend 
on different aspects of the experimental details. The measured resonance cross section 
depends on the energy spread of the ete~ collider, but such finite energy spread does not 
affect the continuum term. On the other hand, the continuum term is sensitive to the 
invariant mass cut in event selection, while such cut hardly affects the observed resonance 
cross section under practical event selection criteria. These are to be discussed in detail 
in the following sections. 


The study of the pure electromagnetic decays of the charmonium states will shed light 
on the understanding of the charmonium decay dynamics. Some theorists, like Suzuki P, 
have tried to probe the interference pattern between one-photon electromagnetic process 
and three-gluon process of the charmonium decays in order to solve the “pm puzzle”. 
For such analysis, the experimental information of those decays of pure electromagnetic 
processes, are of particular importance [L0], since they supply an estimation of the elec- 
tromagnetic part in the decay modes with strong interaction. 


In this paper, we first study the experiment dependence of the cross sections, then 
give an estimation of the cross sections of the pure electromagnetic process without going 
through charmonium resonance based on the form factors, finally the minimum required 
integrated luminosity at off-peak energy point is estimated if meaningful results are ex- 
pected for the decay modes interested. 


2 Experimentally observed cross sections 


2.1 Resonance 


The cross section of the process ete~ — (2.5) — f (where f denotes a certain kind of 
final state) is described by the Breit-Wigner formula 


127 - ely 


opw(W) = W- M22 4 PPM? (2) 


where W is the center-of-mass energy, Te and I’; are the widths of ~(2S) decaying into 
ete” and f, T, and M are the total width and mass of ~(25'). Taking the initial state 
radiative correction into consideration, the cross section becomes [[I]] 


rong ll- 2), (3) 


Orc (W) = f dzF(e, s) 
0 

where s = W?, am = 1 — 8'/s, V's’ is the experimentally required minimum invariant 
mass of the final state f after losing energy due to multi-photon emission; F(x, s) has 
been calculated in many references [I], [2, [3] and H(s(1—~)) is the vacuum polarization 
factor. The radiative correction in the final states are usually not considered [14 [5]. The 
reasons are twofold. In the first place, the hadronic final system is very complicated and 
since the radiative corrections depend upon the details of how the experiment is done, it 
is difficult to give a general, model-independent prescription for them. The second reason 
is that our understanding of the hadronic problem is so crude that there is no need to 
worry about the electromagnetic corrections. In any case, if we find later on that it is 
necessary to do radiative corrections to the hadronic states for some specific problem, we 
can do the calculation then, because the initial state radiative corrections and final state 
radiative corrections can be decoupled to a large extent. 


The ee~ colliders have finite energy spread. The energy spread function G(W, W’) 
is usually a Gaussian distribution : 


GW, W’) = ee a, (4) 


where A is the standard deviation of the Gaussian distribution. It varies with the beam 
energy of the collider. In case of BEPC/BES, A = 1.3 MeV at the C.M. energy of the 
w(2S) [0]. It is much wider than the (2S) intrinsic width of (300 + 25) keV [7]. So the 
experimentally measured resonance cross section is the radiatively corrected Breit-Wigner 
cross section folded with the energy spread function: 


Cenp(W) = J dW'o>..(W'\G(W’, W), (5) 


where Gre, is defined by Eq.(3). 
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Figure 2: Inclusive hadronic cross section of w(2S): ogw for Breit-Wigner cross section; 


Orc. the cross section with radiative correction; Cesp the measured cross section on a 


collider with A = 1.3 MeV. 


In Fig. B, three cross sections are depicted: the Breit-Wigner cross section of Eq.(2); 
the cross section after radiative correction by Eq.(3), and the experimentally measured 
cross section by a ete collider with A = 1.3 MeV. In the calculation of these cross 
sections, the following parameters of ~(25) are used [I] : M = 3.68596 GeV; T, = 300 
keV; Te = 2.19 keV fi, From the three curves in Fig. P|, we see that the radiative correction 
reduces the height of the resonance. It also shifts the peak position to above the w(2S) 
nominal mass. The reduction factor p and the shift of the peak A\/s,__, are approximately 


expressed by 


Max 


ex (1 +6) (6) 
P M , 
AV Smax = i (7) 
where (3 is defined as 
gpa (in a i 
T e 
and ' ; 
3 a(t 1 of 9 T 
p= 3848 (5 5) +8 (5 A (8) 


*Throughout this paper, we use this set of data and A = 1.3 MeV for numerical calculations. 


At the (2S) mass, 6 ~ 0.0779 and 6 = 0.06. So for y(25S), the reduction factor 
p © 0.51 and the shift of the peak AV/s,,,,, © 9 keV. The energy spread further lowers 
down and shifts the experimentally measured 7)(2S) peak. In the case of a collider with 
A = 1.3 MeV, the maximum height of the 7(2S)) peak becomes 640 nb, and the position 
of the peak is shifted by 0.14MeV above the ¥(2S) nominal mass. 


The y*y~ channel deserves special discussion here, since it is parallel to those hadronic 
channels in w(2S) decays which only go through electromagnetic interaction, such as w7° 
and mr. Since this is an exclusive channel, there is interference between the continuum 
and the y(2S) amplitudes. Such interference can be seen clearly from the scan of the 
w(2S), see Fig. B| In Fig. Q, the cross sections of inclusive hadrons and pty~ pairs are 
depicted for comparison. Here in the calculation of radiative correction, the upper limit 
of integration £m in Eq.(3) is taken to be (1 — 4m?/s), which means all allowed phase 
space for multi-photon emission is integrated. The peak of uu curve is shifted more 
than that of the inclusive hadrons, to 0.81 MeV above the W(2S) nominal mass. 


2.2  One-photon continuum 


In the total experimentally measured cross section, the one-photon continuum term has 
different features from the resonance. For the uu channel, this term is expressed in the 


Born order as g 
Ara 


oa (9) 
After radiative correction, this term depends sensitively on the upper limit of the integra- 
tion £m in Eq.(3). At the C.M. energy of the w(2S) mass, the radiatively corrected cross 
section is 8.33 nb, when £m = 1 — 4m, / s. This value means integrating all the allowed 
phase space of multi-photon emission. If £m = 0.2, which means that the final u" u~ 
has the minimum invariant mass of ,/(1 — 0.2) My2s), after losing energy to multi-photon 
emission, then the radiatively corrected cross section is 6.25 nb. In the actual experimen- 
tal situation, the invariant mass cut or its equivalence is usually imposed to remove the 
utp” from J/ decays. For the resonance term, as long as 
T: 


m EFT. 10 
Im > 7 (10) 


O~* Born (s) = 


it is insensitive to £m. Another important feature of the continuum term is that since it 
is a smooth function of s, the finite energy spread of the ete~ collider does not change 
the measured value. 


For the inclusive hadronic final states, the continuum term is expressed in the Born 


order as 7 
Ara 


35 R(s), (11) 
where R(s) is the R-value at C.M. energy ,/s. Since R(s) is a slowly varying function 
of energy, so qualitatively it shares the same feature of this term for pt” final state. 
After radiative correction, it depends sensitively on the upper limit of the integration £m 


O+* Born (s) = 
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Figure 3: Contributions of three parts to the cross section: dashed 
line for QED continuum (of); dotted line for resonance (o®); 
dash dotted line for interference(o‘); solid line for total cross 
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Figure 4: Cross sections in the vicinity of W(2S) for hadron (a) and 


uit ~ (b) final states. The solid line with arrow indicates the peak 
position. 
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in Eq.(3), i.e., on the experimental cuts; but the finite energy spread of the ete” collider 
hardly changes the measured value. 


3 Cross sections of two exclusive electromagnetic pro- 
cesses 


The general discussions in previous section is extended to exclusive processes. For illus- 
trative purpose, we discuss two pure electromagnetic ones, which are similar to ete” > 


wr. 


3.1 Resonance 


In experiments, usually the branching ratio is measured by the events number of certain 
channel (N+) divided by the total number of the resonance events( Nys) ): 

B; = of as Bua a (12) 

Nys) Tys) 

where oy is the measured cross section of the specific channel and oy 25) the total produc- 
tion cross section of (25). With the y (2S) parameters and ete~ collider energy spread 
used in previous section, the total cross section of ~(25) is 640 nb (=oy2s)). The decay 
(2S) — wr’ is reported a branching ratio of (3.6 + 1.6 + 0.6) x 10~° [9]; the other 
decay ~(2S) — ttr" is reported a branching ratio of (8 +5) x 10~° by DASP [Pd], but a 
later result by BES is (8.4 +3.5t1$) x 1076 Bi]. If we check the original references, these 
numbers actually mean that the BES measured cross section of ete~ — wr? at (29) 
mass is 2 x 10~*nb while for ete~ — rtr” it is 5 x 1073nb. 


We notice that the observed cross section of these final states are 3 to 4 orders of 
magnitude smaller than the total inclusive hadron cross section of the continuum process 
which, according to Ref. PJ], is about 15 nb. So it could be that a substantial part of the 
experimentally measured cross section comes from the continuum instead of the w(2S) 
decay. Therefore it is essential to know the production rate of wr? and mtr due to the 
continuum process in order to get the correct branching ratios of Y(2S) decays to these 
modes. 


3.2 One photon continuum 


Here we give an estimation of the possible magnitudes of ete~ — y* — wr? and ete~ = 
y* — ata” by currently available phenomenological models. These two processes are 
calculated by their form factors 


4ra? |p|° 

Cm) = 3 apa Fan(s)’, (13) 
8ra? |p|? 

gal = SEO pa (14) 


where p is the momentum of the final particles, Fi,,(s) and F,,(s) are the form factors 
for wr and m+ respectively. Phenomenological model [IQ] gives 


m? M2, 
For = Pyr — e, 15 
(9) = FO t (15) 
with 
F,,(0) = 2.3 GV. (16) 


At Eom = 3.686 GeV, it gives 3.6 x 10-%nb at Born order. It is almost unchanged by 
radiative correction if £m = 0.2 in Eq.(3)fI 


For ntr, the vector dominance model (VDM) BJ gives : 


m2 


FYPM(S) = —— (17) 


en 
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At W(2S) mass, it gives 3.2 x 10~°nb at Born order, which is almost unchanged after 
radiative correction if £m = 0.9} 


From above estimations, we can see that possibly some of the observed wr? and a 
large fraction of the observed 1*z~ events at (2S) peak may come from the continuum 
process instead of the w(2S) decay. It should be noted that different models for the form 
factors give very different results at the 7(2S) mass region [23], here we only want to get 
a feeling about the magnitude of the continuum contribution. 


For ~(3770), the resonance cross section is only around 8 nb; whereas the hadron cross 
section from the continuum is about 13 nb. To determine its branching ratio of exclusive 
light hadrons, we must know the continuum cross section of these final states. 


3.3 Energy dependence of measurement 


From the foregoing discussion, we know that the behavior of two terms of the cross section 
is rather distinctive with the energy variation: the resonance changes manifestly but the 
continuum part almost imperceptibly. The relative proportion of two terms of the cross 
section, which plays an important role in determination of branching ratio, could be rather 
different at different energy. 


In actual experiment, the beam energy could drift away from the position of the 
maximum inclusive hadron cross section during a considerable long time running, which 
in turn leads to considerable variation of resonance cross section. Taking wtp” final state 
as an example, if the energy drifts upwards or downwards from the (2S) peak position 
of the inclusive hadron by 0.5 MeV, the uu resonance cross section changes from 4.47 
to 4.20 nb, or equivalently the variation is about 6%. 


The radiatively corrected cross section is 3.9 x 107ĉnb if £m = 0.3 and 3.3 x 107° if £m = 0.1. 
tThe radiatively corrected cross section is 3.4 x 107°nb if £m = 0.3 and 3.0 x 107°nb if £m = 0.1. 


Furthermore, for the process whose continuum cross section is comparable with that 
of resonance, its interference cross section accounts for a considerable part in the total 
cross section, and varies violently in the vicinity of the peak position. Also for wt final 
state, its interference cross section changes from —0.14 nb to —0.96 or 0.66 nb respectively 
for the energy drifted by 0.5 MeV lower or higher than the hadron peak position. 


Combining the variation of the resonance and the interference, the maximum change 
of the resonance cross section can reach 1.1 nb, or 25% for utu channel. In other words, 
if we use the result from utu” final state as an approximation, 25% uncertainty should 
be taken into account in the branching ratio estimation for the 0.5 MeV energy driftf}. 


4 Other 7(2S) decay modes 


In previous sections, we have discussed the pure electromagnetic processes in (2:5) de- 
cays. Most other w(2S) decays which strong interaction plays the leading role also have 
continuum amplitude contribution if the ~(2S) is produced in ete™ collision. Although 
the amplitude maybe smaller, but for those suppressed channels, like VP and VT modes, 
the electromagnetic decays and the continuum contribution could play an important role. 
For example, now the upper limit of ~(2S) — pr branching ratio has been pushed down 
to 2.9 x 10-5 B]. Under BES condition, this means that the upper limit of the cross 
section is 2 x 107? nb. The coupling of this channel to virtual photon is one third of w7°. 
Using the estimate of previous sections, its continuum cross section could be 1.2 x 107~3nb 
at (2S). If we push down further the upper limit of this decay, we need to take into 
account the contribution from the continuum. This is also true for the upper limit of 
(2S) = K*+K*(892)- + c.c.. Another VP decay mode, 7)(2S) > K*® K? + c.c., has been 
measured to have a branching ratio of (0.81 + 0.24 + 0.16) x 10~* [B4]. This means that 
BES measured cross section of this channel is 5 x 10~?nb at (29). Theoretically, the 
coupling of this channel to virtual photon is two thirds of w7?, so the estimated continuum 
cross section by the form factors of Eq. (15) is 3.2 x 10~'nb. This has to be considered 
in high precision experiments of the coming generation accelerator and spectrometer, like 
CLEO-c and BES-III. 


In order to know whether the observed suppression of VP and VT modes in 7(2S’) de- 
cays are due to the absence of strong interaction amplitude, or the destructive interference 
between the electromagnetic and the strong amplitudes, or just an incidental destructive 
interference between these two and the continuum process in our particular experiment, 
we need to know their coupling to virtual photon. 


5 Estimation of needed data taking at BES 


As discussed in previous sections, in order to extract the branching ratios of (29) to w7? 
and n*m from experimental data, we need the coupling of these channels to a virtual 


*If the energy drift is about 0.2 MeV, the uncertainty is only 9% correspondingly. 
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photon. This can only be done by the direct measurement of their production cross 
sections nearby but off y(25) resonance. The same is true for ~(3770) decaying into light 
hadrons. To satisfy the need of both w(2S) and ~(3770) physics, the best energy to take 
the data should be below %(25) peak, for example, at Eem = 3.67 GeV. Here two factors 
should be considered: one is the data-taking point should not be too faraway from the 
peak position, so that the cross section of the continuum process is almost equal to that 
at the peak position; another factor is the data-taking point could not be too close to the 
peak position, so that the resonance cross section is small enough to be neglected. Taking 
the inclusive hadron final state as an example, at Eem = 3.67 GeV, the cross section 
variation of the continuum process 
On( Lem = 3.67 GeV) — on( Eem = 3.686 GeV) 


NE 105%, 18 
Sexe on (Ban = 3.67 GeV) a 08) 


At the same time, the ratio of the resonance cross section to that of the continuum process 


Oete-—y(28S)—>hadron 


: ) = x 2.18%. 
continuum / Eem = 3.67 GeV Oete-—y*—>hadron / Pem = 3.67 GeV 


( resonance 


(19) 
That is to say, at Eem = 3.67 GeV, the variation of the continuum process cross section 
and the proportion of the resonance cross section are both fairly small. 


Another reason to take the data below w(2S') peak is that ~(2S) and (3770) are too 
close to find a suitable data-taking point between them. On one hand, w(3770) is a wide 
resonance, the resonance effect can extend far beyond the peak position. On the other 
hand, the effect from the radiative tail of (2S) can reach the vicinity of the (3770) 
peak. For example, at one ['(w(3770)) below (3770) peak, the resonance cross section is 
0.8 nb, which is mostly DD; the radiative tail of (2S) is more than 3 nb; while the cross 
section for continuum process is about 13 nb. That is to say, the radiative tail of ~(2S) 
is too large to be neglected. 


Finally, we discuss quantitatively how much data are needed. 


e The typical electromagnetic process, like w(2S) — w7®, is measured to have a 


branching ratio of 3 x 1075, since the cross section by BEPC/BES at w(2S) peak 
is 640 nb, that means the w7? cross section is 2 x 107? nb. So for every pb~!, we 
obtain 20 produced events. 


If we use the form factor in Eq.(13) to estimate the cross section, at 3.67GeV it is 
3.8 x 107? nb with £m = 0.2 in radiative correction. If so, for every pb~!, we obtain 
3.8 produced events. 


e Another pure electromagnetic process ~)(2S) — ntr" is measured to have a branch- 
ing ratio of 8 x 1076, corresponding cross section is about 5 x 1073 nb. For every 
pb !, we obtain 5 produced events. 


If we use the form factor in Eq.(14) to estimate the cross section, at 3.67GeV it is 
3.3 x 107%nb with £m = 0.2 in radiative correction. If so, for every pb~', we obtain 
3.3 produced events. 
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e ete” — pr is one of the most interesting processes. Each of the three final states, 


pr, ptn, pr gets 1/9 contribution as w7? from pure one-photon process. So 
for every pb~!, we expect 6.7 produced events (which is one of the three final states). 
For continuum process, for every pb~!, we obtain 1.3 produced events. 


From the above estimates, the integrated luminosity of the off resonance data sample 
need to be the same order of magnitude as on 7(2S) resonance. 


According to the experience of BES data taking, the typical integrated luminosity 
at 7(2S') mass is about 8 nb~! per RUN; average 25 RUNs can be taken every day, or 
equivalently 0.2 pb~' per day. As an estimation, using these experiment data and form 
factors, the time needed for events numbers taken at the continuum of different processes 
can be worked out and are given in Table [I] 


Table 1: Events number and time needed 


Events number 1 4 10 25 50 
1.3 5.3 13 33 66 
1.5 6.1 15 38 76 
3.8 15.4 38 96 192 


data taking 
duration 


(day) 


In order to obtain reliable signal, and taking both statistical error and data collection 
time into consideration, we recommend that at least 10 pb~' data at off w(2S) peak to 
be taken. From the above estimates, it gets the statistic to match the new sample about 
20 pb! of Y(2S) taken by BES-II. It will serve the purpose to take into account the 
contribution from one-photon continuum in data analysis. 


6 Summary 


In this paper, we discussed extensively the properties of the observed cross section at res- 
onance in ete” experiment. We studied the resonance and the electromagnetic processes. 
We also pointed out the possible uncertainty due to actual measurement circumstance. 
The electromagnetic final states of w7? and mtr” are from two processes: one is ~(29) 
decay and the other the continuum contribution. Using form factors, we found the contin- 
uum contribution is equivalent to the previously reported cross sections of the resonance 
decay. This indicates that the actual branching ratios of these electromagnetic decays 


might be different from the present reported values. 


In order to obtain the correct results of Y (25) and u(3770) decay, we suggest that at 
least 10 pb~! data be taken at the continuum region (Eem = 3.67 GeV). 
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